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ABSTRACT 

The origin of the diffuse extragalactic gamma-ray background (EGB) has 
been intensively studied but remains unsettled. Current popular source candi- 
dates include unresolved star-forming galaxies, starburst galaxies, and blazars. 

In this paper we calculate the EGB contribution from the interactions of cos- 
mic rays accelerated by Type la supernovae, extending earlier work which only 
included core-collapse supernovae. We consider Type la events in star-forming 
galaxies, but also in quiescent galaxies that lack star formation. In the case of 
star-forming galaxies, consistently including Type la events makes little change 
to the star-forming EGB prediction, so long as both supernova types have the 
same cosmic-ray acceleration efficiencies in star-forming galaxies. Thus our up- 
dated EGB estimate continues to show that star-forming galaxies can represent 
a substantial portion of the signal measured by Fermi. In the case of quiescent 
galaxies, conversely, we find a wide range of possibilities for the EGB contribu- 
tion. The dominant uncertainty we investigated comes from the mass in hot gas 
in these objects, which provides targets for cosmic rays; total gas masses are as 
yet poorly known, particularly at larger radii. Additionally, the EGB estimation 
is very sensitive to the cosmic-ray acceleration efficiency and confinement, espe- 
cially in quiescent galaxies. In the most optimistic allowed scenarios, quiescent 
galaxies can be an important source of the EGB. In this case, star-forming galax- 
ies and quiescent galaxies together will dominate the EGB and leave little room 
for other contributions. If other sources, such as blazars, are found to have impor- 
tant contributions to the EGB, then either the gas mass or cosmic-ray content of 
quiescent galaxies must be significantly lower than in their star-forming counter- 
parts. In any case, improved Fermi EGB measurements will provide important 
constraints on hot gas and cosmic rays in quiescent galaxies. 
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1. Introduction 


The first observation of the diffuse extra, galactic gamma-ray background (EGB) was 
reported by the SAS-2 satellite (Fichtel et al. 1977, 1978). Recently, the Fermi Gamma-ray 
Space Telescope updated the EGB determination from the Energetic Gamma-ray Experi- 
ment Telescope (Sreekumar et al. 1998) and provided the most reliable EGB observations so 
far (Abdo et al. 2009a). Measurements of EGB are model-dependent in that they require 
subtraction of the large foreground emission from our Galaxy (e.g., Hunter et al. 1997). 
The accuracy of the EGB measurement thus greatly depends on our understanding of the 
Galactic emission. Despite the difficulty in its observation, the EGB encodes important 
information about the highest-energy environments in the cosmos. 

The EGB arises from the combination of all the unresolved extragalactic gamma-ray 
sources (e.g., Dermer 2007b; Stecker k Venters 2011). “Guaranteed” EGB components arise 
from unresolved counterparts of known extragalactic populations, namely blazars (those 
AGNs that have their relativistic jets pointing at us; e.g., Padovani et al. 1993; Stecker et al. 
1993; Mukherjee k Chiang 1999; Pavlidou k Venters 2008; Dermer 2007a; Venters 2010; 
Venters k Pavlidou 2011; Inoue k Totani 2009), as well as normal star- forming galaxies and 
starburst, galaxies (e.g., Pavlidou k Fields 2001, 2002; Prodanovic & Fields 2006; Thompson 
et al. 2007; Stecker 2007; Fields et al. 2010; Makiya et al. 2011). Additional EGB con- 
tributions might arise from more exotic sources, such as dark matter annihilation (Silk & 
Srednicki 1984; Rudaz k Stecker 1991), annihilations at the boundaries of cosmic matter 
and antimatter domains (Stecker et al. 1971), massive black holes at redshifts of 2 ~ 100 
(Gnedin k Ostriker 1992), and primordial black hole evaporation (Page k Hawking 1976). 

In this paper we will focus on the EGB contribution from both star-forming galaxies and 
quiescent galaxies. Quiescent galaxies refer to galaxies with little or no active star formation, 
and these objects have not been included in EGB estimations. In terms of galaxy types, 
quiescent galaxies usually include all elliptical galaxies and some SO galaxies. However, the 
important factor for the EGB estimation is not the galaxy type but the amount of star 
formation. Therefore we will separately consider star-forming and quiescent galaxies, and 
assume no star formation in quiescent galaxies. We will not consider the EGB contribution 
from starburst galaxies in this paper, due to the larger uncertainty in the cosmic-ray prop- 
agation in such galaxies (e.g., Thompson et al. 2007: Lacki et al. 2011). Also, recent work 
suggests that starburst galaxies only have small contribution to the EGB (Stecker k Venters 
2011). We follow the criterion adopted in Fields et al. (2010) to distinguish star-forming 
and starburst galaxies. The EGB energy range we consider in this paper is from ~ 30 MeV 
to ~ 30 GeV, which is the energy range covered by Fermi data, and includes the regime in 
which star-forming galaxies may contribute substantially to the EGB. 



Our focus here is on the EGB contribution arising from hadronic cosmic-ray interactions 
with the interstellar medium of their host galaxies, specifically pion production and decay 
pp — > 7T° -> 77. The most favored possibility of the cosmic-ray production sites in galaxies 
are supernovae (SNe). Fermi and air Cerenkov observations detect individual remnants of 
both core-collapse (hereafter CC) (Abdo et al. 2010; Acciari et al. 2010; Weekes et al. 1989) 
and Type la events (Acciari et al. 2011; Acero et al. 2010). The energetics of these objects 
are consistent with the requirements efficient cosmic ray acceleration (e.g., Abdo et al. 2010; 
Reynolds k Ellison 1992), and the GeV spectra of these objects are consistent with pionic 
emission and thus hadronic acceleration (e.g., Abdo et al. 2009b; Tanaka et al. 2011; Ellison 
et al. 2012). Diffuse Galactic emission is similarly consistent with pionic emission dominating 
(Abdo et al. 2009a, c). These data thus give empirical grounding to the long-held belief that 
SNe of all types are the dominant engines that accelerate Galactic cosmic rays (e.g., Baade 
k Zwicky 1934; Ginzburg & Syrovatskii 1964; Ellison et al. 1997). We are thus interested in 
SNe of both types, all of which accelerate cosmic rays in their host galaxies. 

Many groups have studied the EGB emission from cosmic rays accelerated by SNe in 
star-forming galaxies (e.g., Dar k Shaviv 1995; Prodanovic k Fields 2006; Fields et al. 2010; 
Stecker k Venters 2011; Makiya et al. 2011). Some estimations suggest that star-forming 
galaxies can be the dominant source of the EGB (Fields et al. 2010), while other groups 
predict that a major contribution of the EGB comes from blazars (Stecker k Venters 2011; 
Makiya et al. 2011; Inoue k Totani 2009). However, there exist large uncertainties from 
the source inputs. Most of the analyses regarding star-forming galaxies focused on the EGB 
contribution from cosmic rays accelerated by CC SNe and implicitly assume that only these 
events accelerate cosmic rays. We extend the analysis of the EGB from star-forming galaxies 
in Fields et al. (2010) to include Type la SNe as accelerators in the Milky Way and in other 
galaxies. 

CC SNe arise in massive stars with short lifetimes, and thus trace ongoing star formation. 
In contrast, Type la SNe result from thermonuclear runaway of white dwarfs accreting mass 
from their companion stars and hence are related to star formation with some delay time. For 
this reason, observations have shown that Type la SNe exist in both star-forming galaxies 
and quiescent galaxies, while CC SNe are rarely seen in quiescent galaxies (Filippenko 2001; 
Mannucci et al. 2005). Observations have suggested that the intrinsic cosmic CC SN rate is 
about 5 times higher than the intrinsic cosmic la SN rate at redshift 2 < 0.4 (Bazin et al. 
2009). Also, studies suggest that the la rate in a star-forming galaxy is much larger than that 
in a quiescent galaxy; this reflects the distribution of delay times between progenitor birth 
and la explosion, which is weighted towards short delays (Mannucci et al. 2005; Sullivan 
et al. 2006), 
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The efficiency of cosmic-ray acceleration by SNe remains poorly understood but is crucial 
for understanding cosmic-ray acceleration physics as well as SN energy feedback. Theories 
propose that cosmic rays are produced by diffusive shock acceleration in the blast waves from 
SN explosions (e.g., Schlickeiser 1989; Berezhko & Ellison 1999). Current studies suggest 
that ~ 30% of the initial kinetic energy from a SN needs to be transferred to cosmic-ray 
acceleration if we assume that SNe are the dominate sources for cosmic-ray production and 
the nucleosynthesis of lithium, beryllium, and boron in the Milky Way (Fields et al. 2001). 
Also, some theoretical predictions expect the cosmic-ray acceleration efficiency in quiescent 
galaxies is much lower than in star- forming galaxies. Dorfi & Voelk (1996) suggest that only 
< 1% of the total explosion energy goes into cosmic-ray energy in quiescent galaxies. This 
is because a SN blast will have weaker shocks due to the large sound speeds of the hot, 
low-density interstellar medium in an elliptical galaxy. 

Understanding the SN rate and their efficiency in producing cosmic rays is critical for 
studying the EGB contributions from these galaxies. Our observational understanding of 
cosmic SNe will increase significantly when the next generation optical survey telescope, 
the Large Synoptic Survey Telescope (LSST), comes online during the next decade. LSST 
is planning to scan the whole available sky, repeated every ~ 3 days, with unprecedented 
survey sensitivity (Ivezic et al. 2008). The project will observe ~ 10 5 CC SNe per year out 
to redshift z ~ 1 (Lien & Fields 2009) and ~5x 10 2 * 4 Type la events out to redshift z ~ 0.8 
(Bailey et al. 2009). The cosmic SN rate in different galaxy classes can thus be measured 
via direct counting to high redshift with extremely low statistical uncertainty. 

In this paper, we will first describe the general formalism of estimating the EGB from 
cosmic rays accelerated by SNe in both star- forming and quiescent galaxies (§ 2). We will 
then discuss the cosmic Type la rate in each galaxy classification that will be used in our 
EGB analysis (§ 3). The estimations of the EGB contribution from Type la SNe in star- 
forming and quiescent galaxies are presented in § 4 and § 5, respectively. Additionally, we 
discuss the uncertainties in the EGB predictions in § 6. Finally, we summarize the results 
in § 7. 


2. General Formalism 

The formalism we adopt generalizes that of Fields et al. (2010) to account for both 

SN types. Integration of the gamma-rav contributions from each unresolved extragalactic 
source over the line of sight to the cosmic horizon gives the well-known express for the EGB 
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where C y (E e m , z) = dN gamm& / dV com dt em dE em is the comoving luminosity density (or emissiv- 
ity) at rest-frame energy E em , and \dtjdz\ — [(1 +z)H(z)}~ 1 = [{l-\-z)H 0 ^£l m (l + z) 3 + Ha] - 1 
for the standard ACDM cosmology'. We use = 0.27. Qa = 0.73, and H 0 = 70 km s -1 Mpc -1 
from the seven-year WMAP data (Komatsu et al. 2011). 


Because the pionic gamma-ray emission is produced from the interaction between cosmic 
rays and the hydrogen atoms in the ISM of each galaxy, the luminosity density is given by 
the product 

(Lj rtgal) (2) 

of the pionic gamma-ray luminosity L 7 of an individual galaxy, times the galaxy number 
density, appropriately averaged. Our problem then divides into two parts. First, we must 
express a galaxy’s gamma-ray luminosity L 1 in terms of galaxy properties such as SN rate 
and gas content, and relate these to galaxy observables. Then we must construct a luminosity 
function dn ga x/dL 1 for gamma-ray-emitting galaxies. 

We first turn to the pionic gamma-ray luminosity from an individual galaxy. This can 
be written as 


Z/ 7 ( E em ) j F?r c ' !7/ em ) n H dV[ S 
r* 0 -» 77 (Fem) A H 
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where F ff o_ >77 (JS em ) represents a spatial average of the gamma-ray production rate per inter- 
stellar hydrogen atom. The total number A# = f n# dV ism of hydrogen atoms in a galaxy 
is obtained by integrating the number density of hydrogen atom nn over the ISM volume. 
Mu is proportional to the total gas mass M gas in a galaxy and can therefore be expressed as 
Mu — Xu M gas /m p , where Xu is the mass fraction of hydrogen atoms and m p is the proton 
mass. 

We take SNe (of both types) to be the engines of cosmic-ray acceleration; this implies 
that the cosmic-ray flux scales as $ cr oc MscRsn^b- Here RsN.es is an effective SN rate 
weighted by the cosmic-ray acceleration efficiency e, discussed below' . Aesc is the escape path 
length, which quantifies the cosmic-ray confinement in a galaxy. In this paper we assume 
A esc to be universal and constant, which leads to a universal galactic cosmic-ray spectrum 
that is the same as that of the Milky Way. Thus the pionic gamma-ray production rate per 
hydrogen atom in a galaxy should scale as r ff o_ >77 (£’ em ) a # cr oc A esc i?sN.eff- We normalize 
the cosmic-rav spectrum to a known galaxy, which would be the Milky Way in our case, 
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The pionic gamma-ray luminosity of a particular galaxy is thus 
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Because of their short lifespans, the rate of CC SNe (and short-delay la events) traces that 
of star formation: i?sx oc ip. Thus we expect a star-forming galaxy’s gamma-ray luminosity 
to scale as L 1 oc M gas i?$ n oc M sas ip. The new class of Fermi-detected star-forming galaxies 
is consistent with this trend (Lenain & Walter 2011). 


This pionic gamma-ray spectrum always has a peak at E em — m n <>/2, at which the two 
gamma-ray photons inherit the rest-mass energy of the decayed n° (Stecker 1971). At large 
energy, the spectrum shows the same asymptotic index as that of the cosmic-ray spectrum, 
which we take to be 2.75. 


Both CC and Type la events should produce cosmic rays and hence pionic gamma rays. 
Therefore the effective SN rate i?sN,eff hi Eq. (6) is a combination of the effective Type la rate 
■Ria.eff = eia Fi a and the effective CC SN rate Fcc.eff = ecc Fee , where ei a and ecc are the 
cosmic-ray production efficiencies of Type la and CC SNe, respectively. There exist different 
definitions of the cosmic-ray acceleration efficiency in current literature. For example, some 
studies present the efficiency as the fraction of the total cosmic-ray production energy out of 
the total kinetic energy output from a SN (e.g., Dorfi & Voelk 1996; Fields et al. 2001; Helder 
et al. 2010), while other studies define the parameter as the percentage of the energy flux 
that becomes relativistic particles after crossing the shock (e.g., Ellison et al. 2007). Most of 
these definitions describe the fraction of the SN explosion energy transferred to cosmic rays. 
Here, we define the cosmic-ray acceleration efficiency e to be the ratio of the SN baryonic 
explosion energy E SN to the resulting cosmic-ray energy E cr , i.e, e = -E CI /E^- Therefore, if 
we assume all SNe have the same explosion energy and the produced cosmic rays have the 
same energy spectrum, the cosmic-ray acceleration efficiency will be proportional to the total 
cosmic-ray production in a galaxy over the SN rate in that galaxy, i.e., e oc $ cr /(A esc Fsn)- 
For the Milky Way, then, we have R^ eS = C t. mw R\ f w + ccc.mw Rq ™ ■ 

Since we normalized our prediction to the gamma-ray production in the Milky Way 
(Eq. 6). the important factor in the calculation is not the absolute value of e, but the 
difference between the acceleration efficiency e in different SN types (la and CC) and galaxy 
classes (quiescent and star-forming). Specifically, we will need to specify the ratios ej a /ecc 
and cq/cs- Unfortunately, these two fractions are poorly known. Thus for our fiducial 
numerical results, we will take ei a /ecc — 1 and eQ,/es = 1- Furthermore, we are unaware 
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of any evidence for a substantial difference between the cosmic-ray acceleration efficiencies 
between la and CC SNe, and thus we will hereafter drop the la and CC notation in the 
acceleration efficiency to simplify the discussion. However, we will retain the notations of 
the acceleration efficiencies for different galaxy type cq and es in our formalism as an explicit 
reminder that the efficiencies are likely to depend on galaxy environment, as expected by 
some theoretical analyses (Dorfi & Voelk 1996; Hein & Spanier 2008). Further possibilities 
of choosing different cosmic-ray acceleration efficiencies will be discussed in § 6.2. 

Star-forming galaxies contain both Type la and CC SNe. Their pionic gamma-ray 
luminosity density £ 7) s can be calculated by averaging over the galaxy density n ga i axy , 
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In quiescent galaxies, there are almost no star formation. We will assume the star-formation 
rate (and thus the CC SN rate) to be zero in a quiescent galaxy. However, Type la SNe do 
exist in quiescent galaxies because these events occur some time after the star formation. 
Therefore the pionic gamma-ray luminosity density in quiescent galaxies £ 7 ,q only comes 
from Type la events, 
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T?o° (^em) and Rsn are the gamma-ray production rate and Type la event rate in a 
standard quiescent galaxy Q0 for normalization. However, since no gamma-ray emission 
from a (jet-less) quiescent galaxy has ever been measured, we will still adopt the values of 
the Milky Way and estimate the gamma-ray luminosity density for quiescent galaxies as 
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Note that since the gamma-ray production from the Milky Way comes from both Type la 
and CC SNe. (£l em ) needs to be normalized to the total SN rate in the Milky Way 

instead of just the Type la rate. 


The total pionic gamma-ray luminosity density will be a combination of emissions from 
both star- forming and quiescent galaxies, that is, £ 7it0 1 = £ 7i s+£ 7 ,q- The EGB contribution 
from cosmic rays accelerated by CC SNe has been carefully examined in Fields et al. (2010) 
and by other groups (e.g., Stecker & Venters 2011 : Makiya et al. 2011). Here we will focus 
on the EGB contributions related to Type la events. In our calculation, we do not include 
the intergalactic EGB absorption > 30 GeV (Salamon & Stecker 1998). 



3. The Cosmic Type la Supernova Rate in Star-forming and Quiescent 

Galaxies 

Type la SNe do not all trace ongoing star formation, because these events have different 
origins from CC SNe. The prevailing scenarios for Type la SN origin include merging of 
two white dwarfs (double degenerate, Webbink 1984), or a white dwarf accreting from mass- 
overflow of its supergiant companion (single degenerate, Nomoto et al. 1984; Iben & Tutukov 
1984). Both of these scenarios involve white dwarfs merging in a binary system, and thus 
Type la SNe are delayed from the formation of the progenitor stars. For this reason, Type 
la SNe are found in all galaxies, including the quiescent galaxies where there is no longer 
star-forming activity. A complete account of the Type la SN contribution to the EGB must 
therefore include contributions from events in star- forming and quiescent galaxies. 

There are many studies of the comoving cosmic la rate density 1Z\ & — dN\ a / (dV com dt ) 
as a function of redshift, most of which focus on the distribution of delay times (e.g., Scan- 
napieco & Bildsten 2005; Sullivan et al. 2006; Kuznetsova et al. 2008; Dilday et al. 2010b; 
Horiuchi & Beacom 2010; Graur et al. 2011). Our adopted total (all-galaxy) cosmic Type 
la rate is based on current observational data. To find this we use the best-fit redshift and 
delay-time distribution of Horiuchi & Beacom (2010), integrated over all delay times (black 
curve in Fig. 1). Although the number of observed cosmic la events is rapidly increasing, 
the data remain sparse beyond z ~ 1, where the rates are thus poorly constrained. 

The cosmic la rate in different galaxy classes (i.e., star- forming and quiescent galaxies) 
as a function of redshift is still poorly understood. In our calculation, we adopt a constant 
value for the cosmic la rate in quiescent galaxies out to redshift z ~ 2 (red curve in Fig. 1). 
We use a normalization based on observations of la rate per stellar mass in quiescent galaxies 
provided by Sullivan et al. (2006) and a non-evolving stellar-mass function from Pannella 
et al. (2009) (see detailed discussion in § 5). Recently, Type la SNe have been observed in 
galaxy clusters. The observed la rates show little redshift evolution within z < 1 (Gal-Yam 
et al. 2002; Sharon et al. 2007; Graham et al. 2008; Mannucci et al. 2008; Dilday et al. 
2010a; Sharon et al. 2010). Since galaxy clusters are mostly composed of quiescent galaxies, 
these results are consistent with our assumption of a constant la rate in quiescent galaxies. 
Moreover, these measurements show that la rate in clusters is ~ 10~ 13 yr"" 1 Ml 1 , which is 
very similar to the value we adopted from Sullivan et al. (2006) F However, w’hile we expect 
a roughly constant cosmic la rate in quiescent galaxies for moderate redshift, this trend must 
fail at some redshift. To be conservative, we therefore placed an artificial cutoff of the cosmic 


1 We are greatly thankful for the anonymous referee for pointing out the measurements of la rate in galaxy 
clusters. 
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Type la rate at z — 2, beyond which uncertainties in the cosmic Type la rate observations 
remain substantial. 

The cosmic la rate in star-forming galaxies (blue curve in Fig. 1) can thus be obtained 
by subtracting the cosmic Type la rate in quiescent galaxies (red curve in Fig. 1) from the 
total cosmic la rate (black curve in Fig. 1). 

Figure 1 shows the adopted cosmic la SN rate as a function of redshift. Although the 
uncertainty in the rate increases significantly at higher redshift. most of the EGB from Type 
la SNe, like that from core-collapse SNe, arises from events at lower redshift (z < 1) (Ando 
& Pavlidou 2009). For example, in our calculation ~ 50% (~ 70%) of the EGB flux comes 
from sources within z < 1 (z < 1.3). Therefore the choice of the Type la rate at z > 1 
only has a small effect on the final estimation of the EGB. The solid-black curve plots the 
total cosmic la SN rate in both star-forming and quiescent galaxies. The red curve shows 
the cosmic la SN rate in only quiescent galaxies. The blue curve represents the cosmic la 
SN rate in only star-forming galaxies. 

Figure 1 also shows the cosmic CC SN rate (dotted-black curve), which is higher than 
the la rate by a factor ~ 5 at z ~ 0 and increase to a factor of ~ 10 at 2 ~ 1. This 
immediately suggests that we should expect CC events to dominate the star-forming EGB 
signal, with the la contribution at a < 20% level. We will see that this is roughly the case 
for the la contribution from star-forming galaxies, but for la events in quiescent galaxies the 
situation is much more uncertain. 


4. The Extragalactic Gamma-ray Background from Type la Supernovae in 

Star-forming Galaxies 

As described in § 2, the EGB luminosity density £i, a s is dominated by two physics inputs: 
the SN rate in a galaxy, which is associated with the amount of cosmic rays, and the total 
gas mass of that galaxy, which accounts for the total hydrogen targets that interact with the 
cosmic rays. To reflect these two physics inputs, we follow the approach adopted in Fields 
et al. (2010) and rewrite the EGB contribution from Type la events (the first term in Eq. 7) 
as below, 

/"•la _ e S 'fcta.S r MW (p \ \ 

7 ’ S “ e M w R\r + £mw b W em) m p { sas ' s) 


(10) 
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Redshift z 


Fig, 1. — The adopted cosmic SN rate. The solid-black curve plots the total cosmic Type la 
rate; blue curve plots the cosmic la rate in star-forming galaxies; red curve plots the cosmic 
la rate in quiescent galaxies. The cosmic CC SN rate is plotted as dotted-black curve for 
comparison. 
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(12) 


and 7?.i at s = (i?i a ,s %aiaxy.s) is the cosmic Type la rate in star-forming galaxies, as shown in 


Fig. 1. In a star- forming galaxy, the galaxy gas mass M gasS and the galaxy Type la rate 
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i?i a can be related to the star-formation rate in that galaxy, which can be connected to the 
observable Ha luminosity Lna.z of the galaxy by ip(Ln a ,z)/(l M © yr~ x ) = L Ha ,z/(1.26 x 
10 34 W) (Hopkins 2004). Therefore we express the gas mass M gas ,s and the Type la rate R\ & 
in terms of Lh q , z - The corresponding galaxy luminosity function at this wavelength can be 
expressed by the Schechter function (Nakamura et al. 2004). 

At a specific redshift, the gas mass in star-forming galaxies M gaSi s and the star- formation 
rate can be connected by 

Mgas, s = 2.8 x 10 9 Mq (1 + z)- p ( x , (13) 

with f3 = 0.571 and u — 0.714, as shown in Fields et al. (2010). The Type la rate in a galaxy 
can be linked to the star- formation rate via some delay-time distribution A(t), 

/•t(d 

Ru{z) oc / f{t - T ) A (r) dr, (14) 

Jo 

where t(z) is the corresponding cosmic age at redshift 2 . The delay-time distribution A(r) 
gives the probability that a Type la SN explodes a time r after the progenitor’s birth. More 
detailed discussion about the delay-time distribution can be found in Appendix A. The 
galaxy luminosity function at a certain redshift for star-forming galaxies in the Ha band can 
be presented in the form of a Schechter function of 

dn 

dLn a ,z. L*,z \ L*,z 

with a = 1.43 (Nakamura et al. 2004). 

Because Type la. SNe a.re delayed relative to star formation, the Type la rate in a 
galaxy depends on the star formation history of the galaxy via eq. (14). Unfortunately, 
this past star formation rate is not easily determined for distant galaxies even when the Ha 
luminosity LH a ,z is available to give the instantaneous star formation rate. However, we can 
investigate the evolution in two simplified cases: pure luminosity evolution and pure density 
evolution. Pure luminosity evolution assumes that galaxy luminosities evolve with redshift, 
while galaxy density stays unchanged, i.e., L*. z in Eq. (15) has redshift dependence and n*, z 
does not. Pure density evolution assumes that galaxy density evolves with redshift, while 
galaxy luminosity remains constant, i.e., n* jZ in Eq. (15) depends on redshift and L+ z does 
not. The real situation should be bracketed by these possibilities. 

In either limit, we can take advantage of the fact that the well-measured cosmic star 
formation rate is p* — (ipn g& \) ~ (I/H a %ai) ~ This fixes the redshift behavior of the 

product so that in the limit where one factor is constant, the other must have the 

redshift dependence of the cosmic star formation rate. 




4.1. Pure Luminosity Evolution 


In the case of pure luminosity evolution, there is no evolution of the galaxy density. Thus 
in the Schechter function of eq. (15), n*, 2 — n* >0 , and thus all redshift dependence lies in Lua.z- 
Therefore, evolution of the star-formation rate in each galaxy, and hence the evolution of the 
galaxy Ha luminosity Ln a .z, must trace the general evolution of the cosmic star- formation 
rate p*. Under this assumption, we can show that (M gas , s) is independent of the choice of the 
delay-time function (see derivation in Appendix A). When adopting the Schechter function 
for dn d f2 Y ^ one wiU find that ( M gas,s) oc (1 + z)~ fi oc (1 + z)-* with a 

local value of (M gaS) s) 2 =o = 6.8 x 10 9 M© (see Appendix A). 

The predicted EGB from Type la SNe in star-forming galaxies is plotted as the solid 
blue line in the left panel of Fig. 2. For comparison, the dashed blue line shows the EGB 
contribution from CC SNe in star-forming galaxies. The shapes of the dashed blue lines 
trace the results in Fields et al. (2010). However, the normalization of the CC SN curves 
is lower by the fraction of the CC SN rate over the total SN rate (~ 0.8 from Bazin et al., 
2009), which is due to the fact that Fields et al. (2010) have implicitly assumed that CC 
SNe produce all of the gamma-ray emission in galaxies. 

Figure 2 shows that the EGB from Type la SNe is around an order of magnitude lower 
than those from CC SNe, which is due to the lower Type la rate in star-forming galaxies. 
As noted above, this is easily understood as a reflection of the small Ia/CC ratio (Fig. 1). 
The black curve in Fig. 2 presents the total EGB emission from both Type la and CC SNe 
in star-forming galaxies. Note that the total EGB emission from star-forming galaxies turns 
out to be very similar to the prediction in Fields et al. (2010), in which the authors assumed 
that all of the EGB contribution comes from the CC events. The reason is that even though 
we added the EGB contribution from Type la SNe, we also lower the EGB emission from CC 
events estimated in Fields et al. (2010) by the corresponding CC SN fraction. Also, the la 
to CC fraction does not change much within z ~ 1, which is the redshift range where most 
of the EGB signals originate. Hence the EGB contribution from CC SNe is always higher 
than those from la SNe by a similar factor. 

The shape of the EGB curves in Fig. 2 traces the general features of the pionic gamma- 
ray energy spectrum. This is because the observed EGB intensity at a specific energy- 
originated from a combination of sources at different redshifts, as described in Eq. (1). 
Therefore, the redshift evolution of the unresolved sources is smeared out in the energy plot 
and mostly affects the normalization of the EGB intensity but not the spectral shape. 
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Fig. 2.— The EGB from SNe in star-forming galaxies. Results in the left panel assume pure 
luminosity evolution. Results in the right panel assume pure density evolution. The dashed 
blue line shows the contribution from CC SNe: the solid blue line shows the contribution 
from la SNe; and the black line plots the total contribution from both CC and la SNe. The 
Fermi data are obtained from Abdo et al. (2009a). 

4.2. Pure Density Evolution 

For pure density evolution, only the galaxy density evolves with redshift while the galaxy 
luminosity does not. Therefore, the star- formation rate ^ in a galaxy also remains constant, 
and the evolution in the cosmic star-formation rate will purely depend on the growdh of 
the galaxy density. Hence, in the case of pure density evolution, Lhq.z+Az = Thq.z- With 
similar calculations as those in the case of pure luminosity evolution (Appendix A) , one can 
find that (M gaS) g) is also independent of the choice of the delay-time function. Additionally, 
(Mgas.s) cc (1 + zY in the case of pure density evolution. 

Results for the case of pure density evolution are shown in the right panel of Fig. 2. 
Again, the solid blue line and the dashed blue line represent the EGB from Type la and CC, 
respectively. The black line shows the combined gamma-ray contribution from both Type 
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la and CC events. Similar to the results of pure luminosity evolution, the EGB from Type 
la SNe is lower than that from CC events because of the lower Type la rate. Moreover, 
the predicted EGB emission is lower if we assume pure density evolution instead of pure 
luminosity evolution. As discussed in (Fields et al. 2010), this is because a typical galaxy’s 
gamma-ray luminosity L 7 oc ipM gas ; in the pure luminosity evolution case, both factors are 
enhanced at early times, while in the pure density case this nonlinear boost is not present. 


5. The Extragalactic Gamma-ray Background from Type la Supernovae in 

Quiescent Galaxies 


Following a similar procedure to § 4, we will now discuss the EGB from cosmic rays 
accelerated by Type la SNe in quiescent galaxies. We again express the EGB luminosity 
density (Eq. 9) in the following form to describe the physics inputs from the average 
gas mass (M g ^ q) and the cosmic la rate in quiescent galaxies R-i^q, 
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and 1Z U _ Q = (i?ia,Q %aiaxy,Q)- E T nlike the star- forming galaxies, where both (M gas . s ) and 
7Zu,s can be related to the observable Ha luminosity, it is easier to connect both (M gas . q) 
and 7?.i a .Q to the total stellar mass M*,q in a quiescent galaxy. 

For the cosmic Type la rate in quiescent galaxies, we adopt the results of Sullivan et al. 
(2006), which link Type la, rates to M*. q directly. These authors assume a bimodal delay-time 
distribution and decompose the la rate into two groups: the long-delay time and short-delay 
time. In their model, the short-delay time group simply traces the star-formation rate, while 
the long-delay time group has a constant probability for all delay times, i.e., A (r) — constant. 
Therefore the Type la rate in a galaxy can be written as 


i?Ia — A M*,q + B 'll). (19) 

The star formation rate i/j ~ 0 in a quiescent galaxy, thus Ru.q — A M*. q. where A7*.q is 
the total stellar mass created throughout the star-formation history in the quiescent galaxy. 



Sullivan et al. (2006) estimated A = 5.1 x 10 14 vr 1 Mq 1 in quiescent galaxies based on 
measurements of the Type la rate in the Supernova Legacy Survey (SNLS). 

According to the observational results in Pannella et al. (2009), the stellar-mass function 
■ jigr-q of early- type galaxies evolves only slightly with redshift. Therefore we simply assume 
the same stellar-mass function throughout all redshift. Also, we find that the Q shown 
in Pannella et al. (2009) can be roughly fitted by the following function 2 , 


dn 
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where C m — 2.05 x 10 -3 Mpc -3 log -1 M®, fx — 10 10 ' 7 M®, and a m = 0.77. Equation (20) 
(the stellar mass function for quiescent galaxies) and Eq. (19) (Type la rate in a quiescent 
galaxy) give the cosmic la rate in quiescent galaxies 


^■Ia,Q 


dM* t q i?I a ,Q 


dn 

dM+Q 


1.53 x 10 -3 yr -1 Mpc“ 3 * * * 


Here the constant rate follows from the nearly redshift-independent quiescent stellar mass 
function. 

Most of the gas content in quiescent galaxies appears to be in the form of diffuse hot gas 
and can be observed in the X-ray (e.g., Forman et al. 1985; Canizares et al. 1987; Bregman 
et al. 1992). However, large uncertainties exist in estimations of the mass of hot gas. Some 
studies suggest that most of the quiescent galaxies are gas-poor (e.g., David et al. 2006; 
Fukazawa et al. 2006), while other studies imply that there can be significant amount of gas 
in these galaxies (e.g, Jiang & Kochanek 2007; Humphrey et al. 2011). Current theoretical 
models suggest that the gas could extend out to a few hundred kpc from the center. The 
total gas mass of a single galaxy can increase by orders of magnitude if one includes the 
gas at large radii (Humphrey & Buote 2010; Humphrey et al. 2011) Therefore observational 
measures of gas mass can vary even for a single galaxy, depending on -whether or not the 
observations enclose a large enough radius to include gas at large distances. If hot gas 
extends to large radii, this could serve as a reservoir of targets for cosmic rays accelerated 


2 This fitting function is based on observations given in Fig. 7 in Pannella et al. (2009). which contains 

measurements of stellar-mass functions at different redshift bins and environments. We choice the data set 

in the lowest redshift bin of 0.25 < z < 0.55 in medium-dense environment (Logiop ~ —2.75; see Pannella 

et al. (2009) for the definition of environmental density p and more details) to perform a y 2 fitting. The 

reduced-y 2 of our fitting function is ~ 0.22. 
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by Type la. SNe, provided the cosmic rays also reach large radii. This too is uncertain, and 
will depend on the cosmic ray confinement vs escape properties of these galaxies. 

Because of the large uncertainties in the gas content in quiescent galaxies, we illus- 
trate a range of EGB predictions for these objects, based on three different gas amounts. 
Since quiescent galaxies are dominantly early-type, we will use estimates of the gas mass in 
early-type galaxies for the amount of gas in quiescent galaxies. The stellar-mass function 
in Pannella et al. (2009) (Eq. 20) gives the average stellar mass in a early-type galaxies to 
be = 5.17 x 10 u M©, which can be converted to the total gas mass by multiply- 
ing a gas-to-stellar mass ratio M gas /M*, i.e., (M gaSj q) = (M gaS) Q/M* ; Q) The three 

models we adopted for different gas amounts correspond to different gas-to-stellar mass ra- 
tios M gas , Q /M*, Q for early-type galaxies. Table 1 summarizes the gas-to-stellar mass ratios 
Tfgas. q / 1 q for the three models we adopted, as well as the corresponding (M gas> q). Gas 
Model 1 estimates the gas-to-stellar mass ratio based on the stellar-mass fraction of the total 
halo mass from Jiang & Kochanek (2007) 3 * * * * * . This model gives the highest gas-to-stellar mass 
ratio of all three models. Gas Model 2 and 3 are both adopted from David et al. (2006), 
in which the authors reported the gas-to-stellar mass ratio for more luminous (Gas Model 
2) and less luminous (Gas Model 3) early-type galaxies. Note that both of the gas mass 
and the la. rate in a quiescent galaxy are constant with redshift, as a result of assuming 
a non-evolving stellar-mass function, consistent with the early-type galaxy observations in 
Pannella et al. (2009). 

Table 1: Summary of different gas models adopted in the EGB calculations. 


Gas Model 

J Mgas,Q /M*,Q 

(Mgas.Q/ 

Reference 

1 

2.14 

1.11 x 10 12 Mq 

Jiang & Kochanek (2007) 

2 

0.01 

5.17 x 10 9 M q 

David et al. (2006) 

3 

0.001 

5.17 x 10 s M q 

David et al. (2006) 


3 Jiang & Kochanek (2007) found that the average stellar-mass fraction of the total halo mass in early- 

type galaxies is M*/M tot ~ 0.026 or 0.056 based on different, assumptions of the halo mass dynamics. 

Both of these numbers are significantly lower than the cosmological baryon-to-mass ratio Qt,/Q m ~ 0.176 

measured by WMAP (Spergel et al. 2007). If we assume that the baryon-to-mass ratio in a galaxy can 

be well-represented by the cosmological ratio, i.e., M- 0 /M tot ~ Qb/Q m , the result from Jiang & Kochanek 

(2007) implies a large amount of gas mass in early-type galaxies, which can be estimated by M gaSt Q = 

(Mbaryon.Q - M*.q) ~ M*, q ( ~ The values of M*/Af t ot ~ 0.026 and 0.056 correspond to 
-Wgas.Q = 5.77 M*,q and Mgas.Q = 2.14 A7*,q, respectively. Here we adopt the latter number to be more 
conservative in our estimation. 
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The red curves in Fig. 3 plot the EGB estimation from Type la SNe in quiescent galaxies, 
with different line styles correspond to estimations from different gas amounts (solid line: 
Gas Model 1, dashed line: Gas Model 2, dotted line: Gas Model 3). The EGB emissions 
from SNe in star-forming galaxies are plotted as blue curves for comparison. The black 
curves in Fig. 4 plots the total EGB emissions from SNe in both star-forming and quiescent 
galaxies for different gas models. For both Fig. 3 and 4, the left panel plots results under 
the assumption of pure luminosity evolution for the star-forming galaxies. The right panel 
shows the EGB predictions assuming pure density evolution for the quiescent galaxies. 

The estimated EGBs shown in Fig. 3 are linearly proportional to the average galaxy gas 
mass (Mgas.Q) and thus to the adopted gas-to-stellar mass ratio, as required by Eq. (16). Gas 
Model 1 predicts the highest contribution in EGB from quiescent galaxies, where the result 
is significantly higher than the EGB emission from Type la events in star-forming galaxies, 
while Gas Model 2 and Gas Model 3 suggest much lower EGB emission from quiescent 
galaxies. In our estimation, the two important factors that affect the overall EGB are the 
average gas mass and the cosmic SN rate. In general, the cosmic Type la rate is about a 
factor of 5 — 10 smaller than the cosmic core-collapse SN rate. Additionally, the cosmic Type 
la rate in all quiescent galaxies averaging over the entire redshift range is lower than that in 
star- forming galaxies by around a factor of 3. 

As for the gas mass in quiescent galaxies, Gas Model 1 assume the gas amount to be 
about two orders of magnitude higher than that in star-forming galaxies, while Gas Model 
2 and 3 assume gas amount to be similar or one order of magnitude lower than that in star- 
forming galaxies, respectively. Combining these two factors (gas mass and Type la rate), we 
would expect the EGB from quiescent galaxies to be about 30 times larger (for Gas Model 
1), 3 times smaller (for Gas Model 2), or 30 times smaller (for Gas Model 3) than that from 
Type la in star-forming galaxies. Because of the EGB estimation is very sensitive to the 
gas amount, measurement of the EGB could put constraints on the gas mass in quiescent 
galaxies. 


6. The Uncertainties in the Extragalactic Gamma-ray Background Analysis 

6.1. Star-foi'ming Galaxies 

The gamma-ray emission from normal star-forming galaxies is the most constrained 
among the three galaxy classes that are considered possible SN-induced EGB sources (i.e., 
starburst galaxies, normal star- forming galaxies, and quiescent galaxies). The main uncer- 
tainties in the EGB prediction of star-forming galaxies come from four factors, as described 
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Fig. 3. — The EGB from SNe in both star-forming and quiescent galaxies. Results in the 
left panel assume pure luminosity evolution for star-forming galaxies. Results in the right 
panel assume pure density evolution for star-forming galaxies. The dashed blue line shows 
the contribution from CC SNe; the solid blue line shows the contribution from la SNe in 
star-forming galaxies; the red lines shows the contribution from la SNe in quiescent galaxies 
based on different gas models (solid-red: Gas Model 1. dashed-red: Gas Model 2, dotted-red: 
Gas Model 3). The Fermi data are obtained from Abdo et al. (2009a). 

in Fields et al. (2010): (1) uncertainty in the pionic gamma-rav production rate r^^ 77 (£’ e m), 
which is ~ 30% (Abdo et al. 2009c), (2) uncertainty in the normalization of the Galactic SN 
rate , which is ~ 40% (Robitaille & Whitney 2010), (3) uncertainty in the luminosity 
scaling in (M gaSiS ), which is ~ 25% (Fields et al. 2010), and (4) uncertainty in the normal- 
ization of the cosmic SN rate 7 £sn,s- which is ~ 16% resulting from the uncertainties in the 
cosmic CC SN rate TZcc ~ (1-0 ± 0.2) x 10“ 4 yr _1 Mpc~ 3 (Horiuchi et al. 2009) and the 
cosmic la rate 7^i a (0.25 ±0.05) x 10 4 yr 1 Mpc 3 (Horiuchi & Beacom 2010). The total 
uncertainty in the EGB prediction will then be rv 10 ±0 ' 25 . 

The upcoming large synoptic surveys, such as the LSST, will provide novel information 
in both the cosmic SN rate and how they depend on the galaxy types out to high redshift. 
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Fig. 4. — The total EGB from la and CC SNe in both star-forming and quiescent galaxies. 
Results in the left panel assume pure luminosity evolution. Results in the right panel assume 
pure density evolution. The Fermi data are obtained from Abdo et al. (2009a). 

Within one year of observation. LSST is expected to detect ~ 10 5 SNe out to z ~ 1 and 
thus achieve a statistical precision of less than a few percent in the cosmic SN rate (Lien & 
Fields 2009; Bailey et al. 2009). Hence, LSST will essentially remove the uncertainty from 
the cosmic SN rate in the EGB analysis and make the EGB a better tool for studying cosmic 
rays and gamma-ray physics. Moreover, such a large SN population will provide excellent 
statistics for study of how SN rates evolve as a function of SN type (and sub- type), host 
galaxy type and star- formation rate, and cosmic environment. Thus it will be possible to 
directly measure the Type la rate in quiescent galaxies and its evolution with redshift. 


6.2. Quiescent Galaxies: Hot Gas and Cosmic-Ray Propagation 

For quiescent galaxies, many characteristics related to their gamma-ray emissions are 
poorly understood. The Type la SN rate in these systems could be uncertain up to a factor 
~ 2, particularly at z > 1. We have also seen that the gas mass in these systems is even more 




poorly known. Published estimates of the gas content in quiescent galaxies vary by orders of 
magnitude. Observations that use X-rays as tracers of hot gas found significantly less gas in 
quiescent galaxies (e.g., David et al. 2006; Fukazawa et al. 2006). However, several studies 
imply a much higher gas content in quiescent galaxies than previously thought, when using 
gravitational lensing and modeling that involved dark matter (e.g., Jiang & Kochanek 2007; 
Capelo et al. 2010; Humphrey et al. 2011). It is difficult to determine what are the causes of 
the large discrepancy of gas quantity from different analyses, because these studies adopted 
different observational and modeling techniques with different galaxy sample. One possible 
reason for the large variation in gas content measurements might come from whether or not 
one includes gas at larger radii, which are predicted by some theoretical models (Capelo 
et al. 2010; Humphrey & Buote 2010; Humphrey et al. 2011). 

If gas content does extend to larger radii in quiescent galaxies, it is crucial to understand 
the cosmic-ray propagation in quiescent galaxies to determine the confinement volume of 
cosmic rays and how likely they can interact with gas at larger radii. Unfortunately, cosmic- 
ray propagation remains poorly understood. Until now, most of the studies have been 
focused on the Milky Way or spiral galaxies and not so much on quiescent galaxies. For 
spiral galaxies, both observations and theoretical modeling suggest that most of the cosmic 
rays are confined within ~ kpc of the disk (Stecker & Jones 1977; Strong et al. 2000, 2004). 
However, this might not be the case for quiescent galaxies. Hein & Spanier (2008) simulate 
cosmic-ray propagation in elliptical galaxies with presence of diffusion, as well as adiabatic 
losses. They apply their formalism not to SN sources but rather to acceleration due to 
an M87-like relativistic jet; thus their detailed calculations are not applicable to our case. 
Nevertheless, their general finding is that cosmic rays expand into a larger volume in a 
elliptical galaxy than in a spiral galaxy, due to the larger minor axis in a elliptical galaxy. 
Additionally, they argue that adiabatic losses are much more important in elliptical galaxies 
and indeed dominate over escape losses. Hein & Spanier (2008) thus conclude that it is likely 
for elliptical galaxies to be extended gamma-ray sources. 

There are two factors that are important for determining the probability of interactions 
between cosmic rays and gas, especially at large radii where the density of gas and cosmic 
rays are likely to be smaller: (1) the mean free time for pion production for each cosmic 
ray particle r pp-^o = (n gas (Tpp^oc) _1 ~ 1 Gyr (0.1 cm 2 /n gas ), where cr^^o is the cross 
section of pion production and n gas is the number density of gas particles, and (2) time 
r es c that a cosmic ray particle takes to propagate though the galaxy before escape. In the 
limit of cosmic rays travel through a (non-magnetized) galaxy radially, that is, without any 
diffusion, the escape time is very short and probability of interactions between cosmic rays 
and gas particles is small. On the other hand, in the presence of disordered magnetic fields 
extending to large radii, cosmic ray propagation will not be radial but diffusive instead, which 
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will make the escape time much larger and increase the chance of cosmic-rav interactions 
with interstellar medium. 

Clearly there are fundamental uncertainties (and opportunities!) in a realistic treatment 
of elliptical galaxy cosmic rays, but even in our simplistic picture the parameters A esc and 
£q are not well-constrained. Although we treated these two quantities to be the same for 
both Type la and CC SNe in all environments due to limited knowledge, it is possible, and 
even likely, that these numbers are different in quiescent galaxies. In fact, Dorfi & Voelk 
(1996) have suggested that the efficiency in quiescent galaxies is at least 10 times lower than 
that in star-forming galaxies (Dorfi & Voelk 1996), which could lower our prediction of the 
EGB from quiescent galaxies by a factor of 10 or even larger. Tang & Wang (2005) model 
the evolution of SN remnants in low'-density hot media and reach similar conclusion that the 
SN heating in such an environment is subtle and cover a large region because of small Mach 
numbers. Likewise, a smaller escape path length, i.e. , a weaker cosmic-ray confinement, can 
also decrease our EGB estimation in quiescent galaxies. Note that adopting different values 
of e and A esc would change the cosmic-ray spectrum, and hence also change the corresponding 
gamma-rav spectrum. 

To summarize, current studies imply that cosmic-ray acceleration might be more diffi- 
cult in quiescent galaxies than in star-forming galaxies due to weaker shocks in low-density 
and high-temperature environments (Dorfi & Voelk 1996; Tang & Wang 2005). However, 
quiescent galaxies might have larger confinement volume than spiral galaxies because they 
have larger semi-minor radii (Hein k, Spanier 2008). Further observational and theoretical 
study of the gas and cosmic-ray content of quiescent galaxies clearly are needed in order to 
pin down the gamma-ray production of these objects even to within an order of magnitude. 


7. Conclusions 

We have calculated the EGB contribution from Type la SNe in both star-forming and 
quiescent galaxies, extending the work of Fields et al. (2010). For star-forming galaxies, 
most of the gamma-ray emission comes from cosmic rays accelerated by CC SNe. This is 
mainly because there are about 5 times more CC SNe than Type la events in star-forming 
galaxies. We find that the net EGB contribution from both SN types is almost the same as 
the Fields et al. (2010) predictions that only included core-collapse events. Our model allows 
for addition of cosmic Type la explosions, but also includes these in the cosmic-ray /star- 
formation ratio, which we normalize to the Milky Way values. Both factors change by nearly 
the same amount, so that the addition of Type la events causes almost no net change in the 
EGB prediction. 
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We also point out that Type la events in quiescent galaxies make a unique contribution 
to the EGB, because these systems lack CC events. We show that the EGB from Type 
la events in quiescent galaxies is highly sensitive to the gas amount in quiescent galaxies, 
which is still poorly known. Based on different gas models adopted, the EGB from Type 
la SNe in quiescent galaxies can vary from two orders of magnitudes higher to an order of 
magnitude lower than those produced by Type la SNe in star-forming galaxies. Therefore 
quiescent galaxies can be dominant source of the EGB if there exist a large amount of gas, 
as suggested by Jiang & Kochanek (2007) and Humphrey et al. (2011). The measurement of 
the EGB will provide useful constraint on the gas amount in quiescent galaxies. Additional 
uncertainties in the cosmic-ray acceleration efficiency and confinement could also change the 
EGB emission in these systems by a few orders of magnitude. Hence the EGB can also 
provide limits on the these two quantities in quiescent galaxies. 

It is thus important to understand the characteristics of cosmic SNe of all types, in order 
to correctly predict their contribution to the EGB. We conclude that the large SN sample 
provided by LSST will offer critical information about the cosmic SN rate for both CC and 
la events, and their dependence on galaxy types out to high redshift. 

The Fermi detection of the EGB contains crucial information about the extragalactic 
gamma-ray source spectrum. Particularly, it can provide an important probe to the gas 
amount, the cosmic-ray acceleration efficiency and the cosmic-rav confinement in quiescent 
galaxies. With our knowledge about SNe increasing rapidly as future synoptic surveys come 
online, the EGB contribution from SNe in galaxies can possibly be disentangled from other 
source candidates. 
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A. The Delay-Time Distribution of Type la Supernovae and Detailed 
Calculation of (M gas S ) in Equation 10 

The delay time of each Type la SN can differ from ~ 0.1 Gyr to ~ 10 Gyr (Mannucci 
et al. 2005; Scannapieco & Bildsten 2005; Gal-Yam & Maoz 2004; Mannucci et al. 2006; 
Sullivan et al. 2006: Maoz et al. 2011). Observationally, the delay times of Type la SNe are 
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usually studied via comparison between measurements of the cosmic Type la SN rate and 
the cosmic star-formation rate, and are usually described by some functions of delay-time 
distribution, which describes the probability of a Type la event with a specific delay time. 
Current proposed delay-time distributions have included a single power law (e.g., Horiuchi 
& Beacom 2010; Graur et al. 2011), a Gaussian (e.g., Strolger et- al. 2004; Dahlen et al. 
2008), and a bimodal distribution (e.g., Sullivan et al. 2006). The difficulty in determining 
the delay-time distribution mainly comes from the uncertainty in the cosmic la SN rate 
measurements. Fortunately, the value (M gaSi s) is independent of the choice of delay-time 
distribution, as we show in the following derivation. 

Based on the relation between the gas mass M gaSi s and the star-formation rate ip (Eq. 13), 
and also the connection between ip and i.e, ip(Ln a , z)/(l Mq yr -1 ) — LH a ,z/(l-26 x 

10 34 W), the M gaSj s for a galaxy at redshift z can be directly linked to the observable H a 
luminosity L Hq , z by M gas>s = 2.8 x 10 9 M 0 (1 + z) - * 8 ( l2 .,vio» \\ Y- Additionally, the Type 
la rate can also be related to Lh Q jZ by 

Rie.,s oc [ L Ha (t - t) A (r) dr, (Al) 

J o 

where L# Q (t — r) = Z<Ha,z+Az(L T ), which is the H a luminosity measured at some earlier time 
t — r or larger redshift z + Az. Therefore (M gaS) s) can be expressed in terms of Lh^.z, 


t Mgas,s) 


f dLtta,z (1 + A) P {LnaxY (fo L n a ,z+Az (t, t) A(t) dr) 


/ dL]ia,z (fo ^Ha,z+ Az 


(Lt) A (r) dr) 


(A2) 


This equation expresses only the redshift-dependent terms and we will further calculate 
how (M gaS; s) evolves with redshift in the case of pure luminosity evolution and pure density 
evolution respectively in A.l and A. 2. Fields et al. (2010) show's that the local value of 

(Mgas, s)z=o = 6.8 X 10 9 M e . 


A.l. Pure Luminosity Evolution 

In the case of pure luminosity evolution, the star-formation rate in each galaxy traces 
the general evolution of the cosmic star-formation rate. i.e. ~ ■ Therefore 

from the directly proportional relation between the star-formation rate ip in a galaxy and the 
galaxy H a luminosity Z/Ha.z; one can trace the evolution of the H a luminosity via the history 
of cosmic star-formation rate, vdiich is well-known out to redshift z ~ 1 (e.g., Hopkins 2004; 
Hopkins & Beacom 2006, and references therein). That is, -y-— - = 
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Therefore the galaxy luminosity at different redshifts can be found by 

r r P*( Z + Az) _ r [hit - T) 

^Ha,z+Az ^Ha,z • / \ — ^Ha,z . / ,\ 

p*{z) p*(t) 

With this relation, (M gaSt s) in Eq. (A2) can be simplified to 
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Therefore (M gaSi s) is independent of the choice of delay-time distribution. 


(A3) 

(A4) 

(A5) 

(A6) 


The assumption of pure luminosity evolution implies that n* t . in the Schechter function 
remains constant and L* )Z in the Schechter function evolves as p* (Eq. A3). Hence the 
redshift dependence of (M gaS) s) can be further calculated using the Schechter function, 
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(A8) 

« (1 + 2)- 9 (^^0))" (AS) 

where L max is the maximum luminosity for star-forming galaxies, which corresponds to the 
maximum star formation defined in Fields et al. (2010). Galaxies with luminosities greater 
than L max are considered starburst galaxies and are not included in this calculation. Addi- 
tionally, we adopt the cosmic star-formation rate p* described in Horiuchi et al. (2009) based 
on current observations. Note that because the factors related to delay-time distribution 
canceled out, this result turns out to be the same as the one obtained in Fields et al. (2010). 


A. 2. Pure Density Evolution 

In the case of pure density evolution Aho.z+Az = Ln a .z as discussed in § 4.2. Thus in the 
Schechter function, L*, z remains constant while n*. z evolves as p*. With similar calculations 
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shown in the case of pure luminosity evolution (A.l) and adopting the Schechter function for 
s ; we can derive the redshift evolution of (M gas ,s) in the case of pure density evolution: 


(Aigas.s) oc 


f dLna,z(l + z) P (L H a, z)“ (/o Ah q ,z A(t) dr)^^f- 
JdL H *, z (Jo L lia . z A (r) dr) dn g^f 
(l + z)-P (Jj A(t) dr) /dL H q, z (L H q, z )“ +1 
(Jo A(r) rir) / dL Ha , z L Hq , z 
(l + *) -/5 /rfW (W ) w+1 


f rj T T ^-galaxy, S 

J wi>Ha. z ^Ha,z dLu 


oc (1 + z)‘ 


(A10) 

(All) 

(A12) 

(A13) 


Again, because (M gas ,s) is independent of the choice of delay-time distribution, the result is 
identical to the one calculated in Fields et al. (2010). 


REFERENCES 

Abdo, A. A., et al. 2009a, Physical Review Letters. 103, 251101 

— . 2009b, ApJ, 706, LI 

— . 2009c, ApJ, 703, 1249 

— . 2010, ApJ, 710, L92 

Acciari, V. A., et al. 2010, ApJ, 714, 163 

— . 2011, ApJ, 730, L20+ 

Acero, F., et al. 2010, A&A, 516, A62+ 

Ando, S., & Pavlidou, V. 2009, MNRAS. 400, 2122 

Baade, W., & Zwicky, F. 1934, Proceedings of the National Academy of Science, 20, 259 

Bailey, S., Bernstein, J. P., Cinabro, D., Kessler. R., & Kuhlma, S. 2009, ArXiv e-prints: 
0912.0201 

Bazin. G., et al. 2009. A&A, 499, 653 

Berezhko, E. G., & Ellison. D. C. 1999, ApJ, 526, 385 



26 - 


Bregman, J. N., Hogg, D. E., k Roberts, M. S. 1992, ApJ, 387, 484 

Canizares, C. R., Fabbiano. G., k Trinchieri, G. 1987, ApJ, 312, 503 

Capelo, P. R., Natarajan. P., &: Coppi, P. S. 2010, MNRAS, 407, 1148 

Dahlen, T., Strolger, L.-G., k Riess, A. G. 2008, ApJ, 681, 462 

Dar, A., k Shaviv, N. J. 1995, Physical Review Letters, 75, 3052 

David, L. P., Jones, C., Forman, W., Vargas. I. M., & Nulsen, P. 2006, ApJ, 653, 207 

Dermer, C. D. 2007a, ApJ, 659, 958 

Dermer, C. D. 2007b, in American Institute of Physics Conference Series, Vol. 921, The First 
GLAST Symposium, ed. S. Ritz, P. Michelson, k C. A. Meegan, 122-126 

Dilday, B., et al. 2010a, ApJ, 715, 1021 

— . 2010b, ApJ, 713, 1026 

Dorfi, E. A., & Voelk, H. J. 1996, A&A, 307, 715 

Ellison, D. C., Drury, L. O., & Meyer, J.-P. 1997, ApJ, 487, 197 

Ellison, D. C., Patnaude, D. J., Slane, P., Blasi, P., & Gabici, S. 2007, ApJ, 661, 879 

Ellison, D. C., Slane, P., Patnaude, D. J., k Bykov, A. M. 2012, ApJ, 744, 39 

Fichtel, C. E., Hartman, R. C., Kniffen, D. A., Thompson, D. J., Ogelman, H. B., Ozel, 

M. E., k Turner, T. 1977, ApJ, 217, L9 

Fichtel, C. E., Simpson, G. A., k Thompson, D. J. 1978, ApJ, 222, 833 

Fields, B. D., Olive, K. A., Casse, M., k Vangioni-Flam, E. 2001, A&A, 370, 623 

Fields, B. D., Pavlidou, V., k Prodanovic, T. 2010, ApJ, 722. L199 

Filippenko, A. V, 2001, in American Institute of Physics Conference Series, Vol. 565, Young 
Supernova Remnants, ed. S. S. Holt k U. Hwang, 40-58 

Forman, W., Jones, C., k Tucker, W. 1985, ApJ. 293, 102 

Fukazawa, Y., Botoya-Nonesa, J. G., Pu, J.. Ohto, A., k Kawano, N. 2006, Ap.J, 636, 698 
Gal-Yam, A., k Maoz, D. 2004, MNRAS, 347, 942 



-27- 


Gal-Yam, A., Maoz, D., & Sharon, K. 2002, MNRAS, 332, 37 

Ginzburg, V. L., & Syrovatskii, S. I. 1964, The Origin of Cosmic Rays (New York: Macmillan) 

Gnedin, N. I., & Ostriker, J. P. 1992, ApJ, 400, 1 

Graham, M. L., et al. 2008, AJ, 135, 1343 

Graur, O., et al. 2011, MNRAS, 417, 916 

Hein, T., & Spanier, F. 2008, A&A, 481, 1 

Helder, E. A., Kosenko, D., & Vink, J. 2010, ApJ, 719, L140 

Hopkins, A. M. 2004, ApJ, 615, 209 

Hopkins, A. M., & Beacom, J. F. 2006, ApJ, 651, 142 

Horiuchi, S., & Beacom, J. F. 2010, ApJ, 723, 329 

Horiuchi, S., Beacom. J. F., & Dwek, E. 2009, Phys. R,ev. D, 79, 083013 

Humphrey, P. J., & Buote, D. A. 2010, MNRAS, 403, 2143 

Humphrey, P. J., Buote, D. A., Canizares, C. R., Fabian, A. C., & Miller, J. M. 2011, ApJ, 
729, 53 

Hunter, S. D., et al. 1997, ApJ, 481, 205 
Iben, Jr., I., & Tutukov, A. V. 1984, Ap.JS, 54, 335 
Inoue, Y., & Totani, T. 2009, ApJ, 702, 523 
Ivezic, Z., et al. 2008, ArXiv e-prints: 0805.2366 
Jiang, G., & Kochanek, C. S. 2007, ApJ, 671, 1568 
Komatsu. E., et al. 2011, ApJS, 192, 18 
Kuznetsova. N., et al. 2008, ApJ, 673, 981 

Lacki. B. C., Thompson, T. A., Quataert, E.. Loeb, A., & Waxman, E. 2011, ApJ, 734, 107 

Lenain, J.-P., & Walter, R. 2011, A&A, 535, A19 

Lien. A., & Fields, B. D. 2009, J. Cosmology Astropart. Phys., 1. 47 



- 28 


Makiya, R., Totani, T., k Kobayashi, M. A. R. 2011. ApJ, 728, 158 

Mannucci, F., Della Valle, M., k Panagia, N. 2006, MNRAS, 370, 773 

Mannucci, F., Della Valle, M., Panagia, N., Cappellaro, E., Cresci, G., Maiolino, R., Pet- 
rosian, A., k Turatto, M. 2005, A&A, 433, 807 

Mannucci, F., Maoz, D., Sharon, K., Botticella, M. T., Della Valle, M., Gal-Yam, A., k 
Panagia, N. 2008, MNRAS, 383, 1121 

Maoz, D., Mannucci, F., Li, W., Filippenko, A. V., Valle, M. D., k Panagia., N. 2011, 
MNRAS, 412, 1508 

Mukherjee, R., k Chiang, J. 1999, Astroparticle Physics, 11, 213 

Nakamura, O., Fukugita, M., Brinkmann, J., & Schneider, D. P. 2004, AJ, 127, 2511 

Nomoto, K., Thielemann, F.-K., k Yokoi, K. 1984, ApJ, 286, 644 

Padovani, P., Ghisellini, G., Fabian, A. C., & Celotti, A. 1993, MNRAS, 260, L21 

Page, D. N., k Hawking, S. W. 1976, Ap.J, 206, 1 

Pannella, M., et al. 2009, ApJ, 701, 787 

Pavlidou, V., k Fields, B. D. 2001, ApJ, 558, 63 

— . 2002, ApJ, 575, L5 

Pavlidou, V., k Venters, T. M. 2008, ApJ, 673, 114 
Prodanovic, T., k Fields, B. D. 2006, ApJ, 645, L125 
Reynolds, S. P., k Ellison, D. C. 1992, ApJ, 399, L75 
Robitaille, T. P., k Whitney, B. A, 2010. ApJ, 710, Lll 
Rudaz, S., k Stecker, F. W. 1991, ApJ, 368, 406 
Salamon, M. H.. k Stecker, F. W. 1998, ApJ, 493, 547 
Scannapieco, E., k Bildsten, L. 2005, ApJ. 629, L85 
Schlickeiser, R. 1989, ApJ, 336, 243 



- 29 


Sharon, K., Gal- Yam, A., Maoz, D., Filippenko, A. V., k Guhathakurta, P. 2007, ApJ, 660, 
1165 

Sharon, K., et al. 2010, ApJ, 718, 876 

Silk, J., k Srednicki, M. 1984, Physical Review Letters, 53, 624 

Spergel, D. N., et al. 2007, ApJS, 170, 377 

Sreekumar, P., et al. 1998, ApJ, 494, 523 

Stecker, F. W. 1971, NASA Special Publication, 249 

— . 2007, Journal of Physics Conference Series, 60, 215 

Stecker, F. W., k Jones, F. C. 1977, ApJ, 217, 843 

Stecker, F. W., Morgan, D. L., k Bredekamp, .J. 1971, Phys. Rev. Lett., 27, 1469 
Stecker, F. W., Salamon, M. H., k Malkan, M. A. 1993, ApJ, 410, L71 
Stecker, F. W., & Venters, T. M. 2011, ApJ, 736, 40 
Strolger, L.-G., et al. 2004, ApJ, 613, 200 

Strong, A. W., Moskalenko, I. V., k Reimer, O. 2000, ApJ, 537, 763 

— . 2004, ApJ, 613, 962 

Sullivan, M., et al. 2006, ApJ, 648, 868 

Tanaka, T., et al. 2011, ApJ, 740, L51 

Tang, S., k Wang, Q. D. 2005, ApJ, 628, 205 

Thompson, T. A.. Quataert, E., k Waxman, E. 2007, ApJ, 654, 219 

Venters, T. M. 2010, ApJ, 710, 1530 

Venters, T. M., k Pavlidou, V. 2011, ApJ, 737, 80 

Webbink, R. F. 1984, ApJ. 277. 355 

Weekes, T. C., et al. 1989. ApJ. 342, 379 


This preprint was prepared with the AAS LATjrX macros v5.2. 



